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Abstract 

It is demonstrated that in the diffusive superconductor-ferromagnet- 
superconductor (S/F/S) junctions the current-phase relation is practically 
sinusoidal everywhere except in a narrow region near the — n transition. 
In this region the second harmonic dominates the scenario of the — n 
transition. We predict a first order transition for the S/F/S junctions with 
a homogeneous F barrier. However, in real junctions a small modulation 
of the thickness of the barrier may favor the continious — n transition 
and the realisation of the Josephson junction with an arbitrary ground 
state phase difference. 

The performed calculations of the second harmonic amplitude provide 
a natural explanation of the recent contradictory results on the second 
harmonic measurements. 

In usual Josephson junctions ( J Js) at equilibrium the phase difference of the 
superconducting order parameter on the two banks is zero pp. However, the 
situation may be drastically different for JJs with a ferromagnetic intcrlayer 
(S/F/S junctions), where for some intervals of the exchange field h and F layer 
thickness d, the ground state corresponds to the phase difference equal to tt 
("7r-junctions") This phenomenon is related to the damping oscillatory 

behavior of the Cooper pair wave function in a ferromagnet (for more references 
and as reviews see The first experimental evidence of a 0— 7r transition 

in S/F/S (Nb — Cu x Nii_ x — Nb) JJs was obtained by Ryazanov et al. |S] from 
the measurements of the temperature dependence of the critical current. The 
0— 7r transition was signaled by the vanishing of the critical current with the 
temperature decrease. Such a behavior is observed for a F layer thickness d close 
to some critical value d c . In fact, it simply means that the critical thickness 
d c slightly depends on the temperature. The temperature variation serves as a 
fine tuning and permits to study this transition in detail. Recently, thorough 
measurements of the critical current were performed in Nb — Cu x Nii_ x — Nb 
junctions with x = 0.53 and d ~ Tlnm 7 and in similar junctions with smaller 
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x = 0.48 and d ~ Ylnm |Hj . The results were contradictory, since the critical 
current at 0— 7r transition in [7] was zero, while in experiments [S] a small 
critical current was observed. 

In this Letter we elaborate a theory describing how the state is transformed 
into the 7r state. It is demonstrated that the critical current of the S/F/S JJs 
does not vanish at the transition and is determined by the second harmonic 
term in the current-phase relation. This second harmonic contribution decreases 
extremely strongly with the increase of the thickness of the F layer and its 
exchange field. The corresponding estimate for the critical current at the 0— 
7r transition in experiments [7] gives the value well below the experimental 
resolution. On the other hand, for the parameters of S/F/S junctions in Ref. [5], 
the calculated amplitude of the second harmonic is close to the experimentally 
measured value. The 0— 7r transition is discontinuous for junctions with a 
homogeneous ferromagnetic barrier. In real S/F/S junctions the modulation of 
the F layer thickness may provide a contribution to the second harmonic with 
the opposite sign [111,10] • If this mechanism prevails then the 0— 7r transition 
would be continuous. This means that by varying the JJ parameters (e. g. 
temperature) it is possible to obtain the S/F/S junction with arbitrary ground 
state phase difference. We also briefly discuss the thermodynamics of the 0— n 
transition. Note that previously the current-phase relations in S/F/S junction 
were calculated numerically in JT] and analytically for several special types of 
the composite SF-FS junctions and short S/F/S junctions in However 
the theoretical approach to treat analytically the diffusive S/F/S junction with 
the — 7r transition (most relevant to the experiment) was lacking. 

The current-phase relation for JJ is sinusoidal only near the critical temper- 
ature T c [J 

j(<p) = Iisixi(p. (1) 

At low temperature, the higher harmonic terms become more and more impor- 
tant. The calculations of the current-phase relations of the S/F/S junctions in a 
clean limit indeed reveal strongly non-sinusoidal j(tp) dependences at low tem- 
peratures [2], EH, EH- However, in the experiments [JJ, [5] the ferromagnetic 
alloys are used as a F layer, and the dirty limit is more appropriate for the 
description of this case. In such a limit in a normal JJs, if the length d of the 
weak link exceeds the characteristic length £ : of the decay of the Cooper pair 
wave function, the critical current is small j c ~ exp(— d/^) and j(ip) is prac- 
tically sinusoidal jaj- We demonstrate that the second harmonic contribution 
is very small ~ exp(— 2d/£ 1 ). Usually the role of the second harmonic is neg- 
ligible and hardly observable. However, in S/F/S junctions the first harmonic 
vanishes at the — 7r transition and the situation occurs to be very different - 
the contribution of the second harmonic becomes predominant! 
The general current-phase relation 

j(ip) = l\ sin ip + I2 sin 2ip (2) 
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corresponds to the following phase dependent contribution to energy of the JJ 
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If we neglect the second harmonic term, then the state occurs for Zj > 0. Near 
a — 7r transition 7i — » and the second harmonics term becomes important. 
The critical current at the transition j c — \I 2 \ and if li > 0, the minimum 
energy always occurs at tp — or tp — ir (Fig. 1). In the opposite case (I2 < 0) 
the transition from to ir state is continuous and there is region where the 
equilibrium phase difference takes any value < tp < it. The characteristics of 
such a "tp— junction" are very peculiar UJ. 

To describe the properties of the S/F/S junction in the diffusive limit we use 
the Usadel equations Recent studies ^7] revealed a very strong varia- 
tion of j c with the F layer thickness, which implies strong magnetic scattering 
effects ^[j. Assuming the presence of the relatively strong uniaxial magnetic 
anisotropy in F layer we may neglect the magnetic scattering in the plane per- 
pendicular to the anisotropy axes (which mixes the spin up and down Green 
functions) and the Usadel equation for the normal G(x,ui,h) and anomalous 
F (x, ui, h) Green functions in the F layer is (see for example 4 ): 
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where the x axis is perpendicular to the junction plane and the F layer corre- 
sponds to — d/2 < x < d/2, Df is the diffusion constant in the F layer and r s 
is the magnetic scattering time. In the spatially uniform case, the equation J2J 
is equivalent to the one from the Abrikosov-Gorkov theory . The equations 
Q must be completed by the boundary conditions at the S/F interface |2"U] . 
Below we consider two limiting cases: transparent interfaces and large interface 
barriers. Moreover, assuming the normal state conductivity at of the F- layer 
small compared to that of the S-layers, a s >> <jf , we may neglect the influ- 
ence of the F-layer on S-layer, i.e. the Green functions in the left S-layer are 

F s (uj) = Ae lv> / 2 /n , G s (w) = sgn(uj)uj/n, where fi = \juj 2 + |A| 2 (for the right 
S-layer tp — * — tp). 

For transparent interfaces, the boundary conditions 1201 express the conti- 
nuity of the Green functions. At T = T c the equation for F(x, uj, h) is linear 
and may be easily solved 0). Just below T c the nonlinear corrections in are 
small and we may apply the method similar to that, used in the problem of the 
non- linear oscillator |2J. After some calculations we have: 



F(x,ui) = acosh(/cix)+6sinh(A;22 : )- 
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6sinh(3fca;)) . 
(5) 
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where the complex wave- vectors ki and k 2 are determined by the relations 

2 _ , 2 _ fc 2 (5ft 2 ~a 2 ) _ (3a 2 +b 2 ) , 2 _ , 2 _ fc 2 (5a 2 -b 2 ) _ (3h 2 +a 2 ) , 
ft l K - 8 4-D /Ts ' ft 2 re - 8 4Dj.r s ' ±ue 

non-linear effects make k\ and ki different from the wave- vector k 2 (uj,h) — 
2(ui+ih+sgn(uj)/T s ) / D f in the solution of the linear equation. The coefficients 
a and b are determined from the boundaries conditions and the current-phase 
relation may be then directly obtained from the formula for the supercurrent 

In the limit h » T c and in the absence of magnetic scattering the Cooper 
pair wave function in ferromagnet decays and oscillates at the same character- 
istic length £f = \fhjD~f [221 • The magnetic scattering leads to the decrease 

of the decaying length £ 1 = £f/y/a + ^l + a? and to an increase of the os- 
cillating length f 2 = £// V Vl + a 2 — a, where the dimensionless parameter 
a = l/QiTg). Note that the product £^2 = £/ is a constant and independent 
on a. If the F layer thickness d > ^ the expression for the first harmonic reads 



.2 /1 1 „.2\V4 
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where j = 47rX , c eSW(0)I)//£ / , the angle 7r/4 < * < tt/2 is determined by 
the relation tan ^ = a + \/i + a 2 , and 5 is the area of the cross section of the 
junction and iV(0) is the electron density of state per one spin projection. The 
amplitude of the first harmonic © reveals an oscillatory decay as a function of 
the F layer thickness. The sign change of 1\ signals the - 7r transition. At 
T T c the critical thickness separating and n phases is d™(T c ) = (jm — VP) £ 2 , 
which decrease with the decreases with temperature 



The second harmonic term is much smaller 

(d/Sf) sin (2d/£ 2 



~Jo^exp(-2d/e i ; 



5 sin (2d/£ 2 + *) + 3 sin (2d/£ 2 - 3*) 



4^/2 8^(2*) 

(8) 

and compared with I\ it contains, in addition to the factor A 2 /T c 2 , an ex- 
ponentially small term exp(— d/^). The amplitude of hid) also reveals the 
oscillatory- like dependence similar to I\(d) but with the decaying and oscil- 
lating lengths two time smaller. It may be directly verified that at the - 7r 
transition (i.e. when d = d") I2 is always positive. 

Now let us consider the limit of large S/F interface barriers which is charac- 
terized by the parameter j B , related to the S/F boundary resistance per unit 
area by 7 S = Rb<Jf. In this case the boundary condition reads j B (dF/dx) x=d , 2 
= F s (x = d/2) G 2 (x = d/2) [201 - Performing the similar analysis as in the per- 
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feet transparency case (7^ — 0) we obtain in the limit d > ^ : 

^=jo(4r) ~, mf. 1/2 eM~d/^)sm(^-d/^). (9) 
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The second harmonic term at a ~ 1 is of the order I2 ~ jo f^"J ex P(~ 

and also oscillates (with a period 7r£ 2 , which is two times smaller than that of I\ ) 

and similar to the transparent interface case is positive at the - 7r transition. 

For small a « 1, the first - tt transition occurs at the F layer thickness 
smaller than £^ [22j • -"- n the considered case (assuming h,T~ l » T ) the first 
- 7r transition occurs at d = d? c = £ ^V3a and 

/ Z f \ 2 A 2 d° c -d 

The second harmonic term at the - n transition is also positive and I 2 ~ 

jo T%/a • ^ ne f° rmma (|9I10I) are written for T <T C but the correspond- 

ing analysis is easily generalized for all temperatures. Besides the change of 
numerical coefficients the expressions for I\ and I2 remain the same. 

Now we demonstrate how the obtained results permit to understand the 
controversy in the experimental search of the second harmonic 0, [SI- From the 
thickness dependence of the critical current in the series of Nb — Cuo ^Ni 53 — 
Nb junctions J7j,JSI we ma y estimate ^ w 1.4 nm and £ 2 ~ 4.1 nm. This gives 
Zj « 2.4 nm, the magnetic scattering parameter a « 1.3 and the exchange field 
h « 600 K. In the experiments [7] the current-phase relation was measured near 
the second - 7r transition at d « 22 nm (the first transition occurs at d « 11 nm 
[T%]). Therefore we may roughly estimate /2/ii ~ 0.1 exp(— ~ 10~ 8 , which 
gives very small value for I2 < 10~ n A, well below the experimental threshold. 
On the other hand the corresponding estimate for the first - tt transition at 
d ~ 11 nm is much more favorable for i 2 observation: I2/I1 ~ 10 -4 , and I 2 ~ 
10~ 6 A. Therefore it would be interesting to perform the similar measurements 
on the junctions revealing temperature mediated first - 7r transition. 

In similar junctions, but with smaller Ni concentration x — 0.48 the second 
harmonic at the - 7r transition was reported for the F-layer thickness d « 17 
nm at 1.1 K (Hj. This is the first (as a function of d) - tt transition. For 
junctions with x — 0.48 ^B] we may roughly estimate £j « 4 nm , £ 2 w 9 
nm, the magnetic scattering parameter a ~ 0.9, nm and ft, 100 if. 

Extrapolating the expressions JfJ) , {HI for low temperature we have for the ratio 
h/h ~ 0.1 exp(— d/^) ~ 10~ 3 which is close to the observed value 3 • 10 _3 [5]. 
On the other hand, for the F-layer thickness d w 19 nm the - tt transition 
occurs at the temperature 5.3 if and the second harmonic was too small to be 
observable. Smaller A/T ratio and larger F layer thickness gives this case less 
favorable for the second harmonic observation. 

We have demonstrated the presence of a small intrinsic second harmonic at 
the —tt transition in S/F/S junction with a uniform barrier. However, there 
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is another mechanism of the negative second harmonic generation due to the 
inhomogeneity of the F layer thickness [5] . Indeed the roughness of F layer in 
the real S/F/S junctions [7], JUiE! is of the order of 1 nm. This means that if 
the characteristic length 51 of the thickness variation (along the contact surface) 
is larger than d, the critical current will vary locally too. On the other hand, if 
61 is much smaller than the Josephson length Aj, which for the current density 
10 6 A/m 2 JTJ is of the order of the junction dimension (50x 50 /iro 2 in J7J)> the 
measured characteristics of the junction will be effectively averaged. At the — tt 
transition we deal with a system where the local current density is alternating 
±Ji and I\ — 0. The resulting local phase variation leads to the appearance 
of the negative second harmonic in the averaged current-phase relation |10|.|5| 
I 2 ^ — \Ii\ (51 / \j) 2 , where Aj is the Josephson length corresponding to the 
current density I\. The 1 nm roughness of the F layer in the experiments 
[5].|16| permits to estimate for d « 17 nm the value |ii| ~ 5 • 10 6 A/m 2 and 
A,/ ~ (10 — 100)fim. At the present time there is no information about the 
characteristic length 51 of thickness variation in the studied S/F/S junctions. 
Taking it as 1 am for the lOx 10 am 2 junction 5H]>!B] we have h ~ — 5(10 2 — 10 4 ) 
A/m 2 while the experimentally observed value is ~ 3 • 10 4 A/m 2 and the sign 
of ]2 is unknown. 

Let us now briefly discuss the thermodynamics of the —tt transition. Near 
the transition temperature T„, the amplitude of the first harmonic Ii may be 
considered as a linear function of T, i.e. I\ — n(T — T n ), while the second har- 
monic term as a temperature independent. The phase dependent contribution to 
the free energy of junction © being Ej(ip) — \—fi (T — T n ) cos ip — -j- cos 2(p] .For 
h > the transition occurs to be I order and at T > (0 phase) 5F = 
-^( T - T -)' while at T < T w , in the tt phase 5F* = -^/x(T OT -T). 
Therefore the latent heat of the transition is q = ^fiT^. Taking the parame- 
ters of the S/F/S junctions [U,EI], we ma ^ es timate /xT^ ~ 10 -4 A and then 
q ~ 3 • 10 -20 J. The S/F/S junctions studied in ^Hj with the - tt transition at 
d « 11 nm reveal the parameter /xT^ and consequently latent heat, which must 
be two order of magnitude larger. 

In the case I2 < the - tt transition is continuous and in the interval 
— |72| < I\ < 1-^2 1 the equilibrium phase difference is determined by cos^ = 
fi(T -T w ) /\I 2 \. The specific heat of this ip- phase is 5C = §^T^. Therefore 
we may expect on experiment an increase of the specific heat by 5C in the narrow 
temperature region T v — fi I/2I < T < + u \I 2 \. If we suppose that the - tt 
transition observed in |S] is continuous , we may estimate 5C in the (p— phase. 
Taking the experimental value of I2 ~ 3 a A we have 5C ~ 10 -18 J/K = aJ/K. 
The recent precise measurements of the specific heat of the superconducting 
microrings demonstrated the possibility to register a specific heat variation 
of the order of 0.1 aJ/K per ring. 

It may be easily verified that (independently on the sign of I2 and then 
on the scenario of the transition) the critical current varies near the transition 

temperature (|7i| << |J 2 |) as I C (T) = \I 2 \ (l + ^ T - J^ \ . The temperature 
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dependence is linear and the slope dI c /dT near T n is l/v2 time smaller than 
that far way from T„. 

At T = T v the current-phase relation is j (ip) — ^sa^^ip and then the 
AC Josephson effect would imply the frequencies two time larger u> — » 2lo — 
2 (^ V) . Note that this circumstance was responsible for the observation of the 
half-integer Shapiro steps at - n transition in experiments [5j and provided an 
additional proof of the nonzero critical current at T — T^. The SQUID with 
such junctions would have the periodicity $o/2 on the magnetic flux. 

The structure of the soliton (Josephson vortex) is rather peculiar for the 
long (along y— axis) junction near a - ir transition. For I2 > 0, the phase 
distribution is determined by the following equation 

d 2 ip 1 . . 

—2=-2-{esmcp + sm2ip), (11) 

where A^o = c^qS/ (8ir 2 t \h\), £ = h/h, and t is the effective junction thick- 
ness. The soliton-type solution is 

( 2(l + e) \ 

tp = arccos 1 ?, f , r • 12 

V l + ecosh 2 (VTTi2//A. /0 )/ 

The variation of the shape of the soliton is presented in Fig. 2. Approaching 
the transition the central part of the soliton with ip rs it grows and finally at 
Ii = the system has two degenerate ground states ip — and ip = tt. If I2 < 0, 
the ip— junction is realized in the interval — I/2I < Ii < \h \- The S/F junctions 
near the temperature T — could provide an excellent possibility to study the 
unusual properties [UliEl of these junctions. 

In summary, we present an analytical solution of the problem of the second 
harmonic contribution to the current-phase relation of the S/F/S JJs in diffusive 
limit in the presence of uniaxial magnetic scattering. Note that very recently 
the case of the isotropic magnetic scattering have been studied numerically in 
|24| with qualitatively similar results. An important conclusion of our work is 
that the - 7r transition is discontinuous for the S/F/S JJs with a homogeneous 
F barrier but may be continuous in real junctions with modulated F layer thick- 
ness. In the latter case a very special ip— junction exists in the transition region. 
The modern microcalorimetric technique could be used for the experimental 
study of the thermodynamics of the - 7r transition and determines its type. 
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Figure captions 

Figure 1. Schematic plot of the phase dependent JJ's energy. The case 
li > corresponds to the discontinuous — ir transition while for I 2 < 
the minimum energy is reached at < ip < ir. 

Figure 2. The change of the form of the Josephson vortex at the discontin- 
uous — 7T transition. The parameter e — I\ / h vanishes at the transition. 



9 



